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Abstract 



We show that a minimal model for viscous fingering with a nematic liquid 
crystal in which anisotropy is considered to enter through two different viscosi- 
ties in two perpendicular directions can be mapped to a two-fold anisotropy 
in the surface tension. We numerically integrate the dynamics of the resulting 
problem with the phase-field approach to find and characterize a transition 
between tip-splitting and side-branching as a function of both anisotropy and 
dimensionless surface tension. This anisotropy dependence could explain the 
experimentally observed (reentrant) transition as temperature and applied 
pressure are varied. Our observations are also consistent with previous exper- 
imental evidence in viscous fingering within an etched cell and simulations of 
solidification. 

PACS number(s): 47.54.+r, 47.20.Hw, 61.30.-v, 47.20.Ky 
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I. INTRODUCTION 



Interfacial instabilities arise in a wide variety of contexts, often of applied interest, such 
as dendritic growth, directional solidification, flows in porous media, flame propagation, 
electrodeposition, or bacterial growth 0. Notwithstanding this disparity, there has been 
a search for unifying common features concerning the more fundamental problem of their 
underlying nonequilibrium dynamics. One such feature seems to be the role of anisotropy 
in determining the observed morphology. 

Thus, the finding that anisotropy is necessary for the needle crystal to solve the steady 
state solidification problem (see e.g. 0) and indeed a critical amount of it to stabilize its tip 
and (possibly) generate side-branches in related local models || motivated the inclusion of 
anisotropy in viscous fingering experiments, either by engraving the plates M, or by using a 
liquid crystal In turn, these experiments in anisotropic viscous fingering confirmed the 

existence of a tip-splitting / side-branching transition controlled by anisotropy and driving 
force. Again, theoretical work on the solidification problem || and numerical integration 
of its (nonlocal) dynamics [10| showed the anisotropy in the surface tension to control the 
transition in this system together with the dimensionless undercooling. To our knowledge, 
no analytical or numerical work on anisotropic viscous fingering [|TTJ has focused on this 
transition so far. 

In this way, a picture that some kind of anisotropy can control the tip-splitting or side- 
branching behavior of different systems emerged. However, it is still not clear how different 
these anisotropics and systems can be. On the one hand, not all kinds of anisotropy seem to 
control the transition. The channel walls in a viscous fingering experiment, for instance, are 
known to play the same role as surface tension anisotropy in free dendritic growth as far as the 
existence of a single finger steady state solution with surface tension is concerned. Moreover, 
even with isotropic surface tension, this steady finger is stable up to a certain critical amount 
of noise but no side-branching is observed, in contrast to free dendritic growth. Above 



this threshold the tip splits, but again no side-branching is observed, in spite of the anisotropy 
due to the channel walls. It is necessary to introduce some other type of anisotropy to 
observe the transition to side-branching. On the other hand, some types of anisotropy not 
directly acting on the surface tension jH| seem to actually control the transition in some 
systems, as is the case of liquid crystal viscous fingering experiments, which varied mainly 
the anisotropy in the viscosity, as well as etched cells, where the exact effect of the grooves 
on the free boundary equations is unclear. Therefore a connection between surface tension 
anisotropy (seen to control the transition in the solidification problem both experimentally 
and in simulations) and other types of anisotropy clearly lacks. Here we present such a 
connection for the case of a simple model of a liquid crystal. 

Specifically, we show that two different viscosities in two perpendicular directions (in 
addition to some already anisotropic surface tension) can be mapped to a two-fold surface 
tension anisotropy (times the rescaled original anisotropic surface tension) through a conve- 
nient axis rescaling. Moreover, we integrate the resulting problem to confirm the existence 
of the morphological transition also for the viscous fingering equations. The numerics use a 



previously developed [14] and thoroughly tested phase-field model for viscous fingering. 

We do find such a transition as a function of the amount of anisotropy and the value 
of the dimensionless surface tension itself (i.e., the driving force for fixed surface tension). 
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The results are consistent with experimental results in viscous fingering with a liquid crystal 
f^H and etched cells [|J , and also with theory and simulations for the solidification problem 

0JEJ- 

The layout of the rest of the paper is as follows: in Sec. |I| we refer to the special features 
of liquid crystals concerning viscous fingering and present a simple model for them. We 
then map this model onto the basic Saffman- Taylor problem with a two-fold anisotropy in 
the surface tension. In Sec. |TJ we briefly describe the phase-field model used and present 
the numerical results. Finally, in Sec. |V| we discuss their consequences for viscous fingering 
with a liquid crystal and consistency with related problems. 



II. MODEL 

In the nematic phase of a liquid crystal its molecules are locally oriented, giving rise 
to anisotropy in the viscosity and surface tension. The degree of orientation depends on 
the proximity of the other phase(s), namely the isotropic (and for some liquid crystals the 
smectic), i.e., it still depends on temperature, and so does the anisotropy, mainly in the 
viscosity |7j. Therefore one should be able to explain the tip-splitting / side-branching 
transition as a function of temperature in the nematics by means of the anisotropy in the 
viscosity alone. 

In a viscous fingering experiment, the director forms a small angle with the velocity 
field, except maybe for the neighborhood of the interface, where it might follow its normal 
direction. So, as a first approximation, one can consider that there is flow alignment, and, 
therefore, a velocity dependent viscosity, which would make the flow nonlaplacian. However, 
in the vicinity of a finger tip we can approximate the direction of the flow for that of the 
finger, so that we can make a minimal model with only two different viscosities: one in the 
direction parallel to the finger and one in the perpendicular one. More details can be found 
in Ref. §. 

Let us now review the formulation of the Saffman-Taylor equations to account for those 
two different viscosities in two perpendicular directions x and y. We will do it for the chan- 
nel geometry, although the result also applies to the circular cell used in the experiments of 
Ref. 0] with minor changes, (i.e. two different viscosities in the radial and tangential direc- 
tions would also map to the standard viscous fingering equations and the same functional 
dependence for the surface tension anisotropy). For the sake of generality, we consider both 
the displaced (1) and the injected (2) fluid to have a certain distinct viscosity (/ii,/^)- 

As in the usual Saffman-Taylor problem, in each bulk we assume the flow to be incom- 
pressible, 

V-m = (2.1) 

(where u is the fluid velocity in the reference frame moving with the mean interface at V^, 
the injection velocity), and also Darcy's law to hold, but now for two different viscosities 

u x = — — d x p (2.2) 

u y = — — (d y p + Pig e ff) - Voo, 
Vy,i 
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where i = l,2 stand for each fluid, the x, y components of u, p is the pressure, 

r] Xji = (12/6 2 )/i xi is an inverse mobility in the x direction, r\ y ^ = (12/b 2 )p, yti , in the y 
direction, p i} the density, and g e ff, the effective gravity in the plane of the channel. Also 
as in the usual Saffman-Taylor problem, on the interface the normal velocity is continuous 
and equals that of the interface, 



r • U\ — r ■ «2 = v n 



(2.3) 



(where r is a coordinate perpendicular to it increasing towards fluid 1 and v n its normal 
velocity), and the pressure has a jump given by Laplace' law, 



Pl -p 2 = <t(</>)k, 



(2.4) 



with a(4>) the (anisotropic) surface tension and k the interface curvature. 

Due to Eqs. (|2.1| ), and ( j2.3|) , the flow can be described by a scalar field ip, the stream 
function, defined even on the interface by u x = d y ip, u y = —d x ip (see e.g. Refs. [[TJ],|I(J). 
However, because of the different viscosities in the x and y directions, the problem is non- 
laplacian (there is vorticity) in the bulk: 



V 2 ^ 



|V x u\ ^ 0. 



(2.5) 



To circumvent this, we rescale the x and y axis by a different factor. We also take 

so 



advantage to adimensionalize the resulting equations in the same way as in Refs. |14| , |T6 



that they can be compared to those in these references, and especially to Ref. |14[ in order 
to generalize the phase-field model described there to the case of anisotropic viscosity. Thus, 
we perform the following change of variables: 



V = a y y, 
W~ 



(2.6) 



where tildes denote new variables, a x , a y have units of length, is a velocity, and W, the 
channel width. We find 



(2.7) 



so that the flow is still incompressible and we can define a new stream function xjj = 
{W/U*) [ip/(a x a y )], which will be laplacian in the bulk if and only if [see Eq. (|2.5|) 1 the 
velocity field is potential in each fluid, 



u 



W 



1 



Vp + a y pigy 



-y 



which is now the case as long as we choose a x , a y to be such that 



(2.8) 



a 2 x Vx,i = CL 2 y rj y ^ = W 2 f]i. 



(2.9) 
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On the interface, Eq. (2J3) will be formally unchanged as long as the choice of a x , a y is the 
same at both sides. Note that, according to Eq. (|2.9|) , this implies that the ratio m = rj x /rj y 
must be the same for both fluids. In an air-liquid crystal experiment, this is obviously not 
the case, but, in the limit in which the viscosity of the air is negligible compared to that 
of the liquid crystal, the anisotropic character of the air viscosity in our model becomes 
irrelevant. In terms of the stream function, Eq. (|2.3|) for the new variables then reads 



dsipi = d s ip 2 



(2.10) 



where s is the arclength along the interface and such that s x f = x x y. 

As for drip, the boundary condition for it will be given by that for u s = s ■ u. Indeed, it 
will have a jump on the interface due to the fact that u is not potential on the very interface 
[see Eq. (|2.8|| because of the jump in fji, which gives rise to a singular vorticity on it: 

(jk + V2)(u- sA - M g2 ) + (fji - 7fr)(us,i + U i>2 ) 



w 



w 2 



dsipi - V2) 



dy 



P2, 



{fji - rf 2 )V oc 



y ■ s 



(2.11) 



and therefore, making use of Eq. 



dfipi - dfip2 = us i - u s ,2 



2 



W 2 (7] 1 + 7] 2 

-c (dfipi + d r ip2 



— d~ s [a{<P)WK] + 



a y gijh - Pi) W_ cV 

W (?7l + 772) Oy 



y ■ s 



(2.12) 



where c = (fji — 772) / (?7i + Now choosing a y /a 2 x = 1/W, Eq. ( p.9| ) yields m = a y /W, 
and defining U* = cV^/m + [mg(p\ — p 2 )] / (Vi + V2) we recover the usual result for viscous 
fingering in a channel (see Refs. fl4Ul6| ), 



dfipi - dfip 2 = -2ds[B(cp)Wn} —2y-s — c (d?ipi + d?ip< 



(2.13) 



with B((p) = cr(<fi)/ [W 2 (fji + 772) Z7*] , except for the m factors in the definition of Z7* — and 
therefore in B(<p) — and the clue fact that Wk and <r(0) are still in the old variables and 
must be rescaled: 



K 



dx 2 



1 + 



dx 



ay_d_y_ 

al dx 2 



1+ 



a y dy^ ~ 

cXi y dec 



-3/2 



-3/2 



l_d^y_ 

Wdx 2 



dy\ 2 



1 + m Ui 



-3/2 



(2.14) 



so that we obtain 



Wk 



1 + (dy/dx) 2 3//2 
1 + m (dy/dx) 2 _ 



K 



1 + (m — 1) cos 2 



3/2 • 



(2.15) 
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where <fi is the angle from x to f. 

To summarize, we recover the usual viscous fingering equations, including Eq. ( [2.13 ) 
which finally reads 



dfipi — dfi>2 = — 2d s [B((J))K} —2y-s — c (df^i + df^2 
but now with an anisotropic dimensionless surface tension of the form 



B(4>) = B x E(0) x 



1 



1 + (m — 1) cos 2 



3/2 



(2.16) 



(2-17) 



where B is the dimensionless surface tension of isotropic viscous fingering 

0o 



Bn = 



W 2 [(fjt - r] 2 )V oc /m + mg(p 1 - p 2 )\ 



(2-18) 



except for the m factors, with cr(0) = <ToS(0). This means that, even with an originally 
isotropic surface tension, that in the rescaled problem has a two-fold anisotropy with a very 
specific form given by the last (third) factor at the r.h.s. of Eq. ( |2.17 ). On the other hand, 
the possible original anisotropy in the surface tension will change its functional form in the 
rescaled problem according to 



E(0) = E(0) = E 



arctanf v^tan 



(2.19) 



and the rescaled problem will have the two-fold anisotropy of the mentioned last factor 
superimposed to the transformed anisotropy of Eq. (|2.19|) [second factor at the r.h.s. of 
Eq. (|2.17|) 1- A similar result was found in a different context, namely for the nematic- 
smectic B transition, where two different heat diffusivities in two perpendicular directions 
could be mapped to the same type of anisotropy in the surface tension and the same type 
of transformation in the original anisotropy [|T7|]. However, note that here the assumption 
is that the growth is in the direction of lowest viscosity (because of flow alignment of the 
director), which results in growth in the direction of largest surface tension (0 = tt/2), 
whereas in Ref. |17j the situation is just the opposite: growth was found to be in the 
direction of lowest diffusivity because that is the direction of lowest capillary length. Also 
for isotropic diffusivities it is known that steady needle crystals can only grow in the direction 
of minimal capillary length ||, although there the anisotropy is assumed to be four- fold. 

Finally, for the described minimal model the original anisotropy in the surface tension 
would be two-fold, e.g. 



S(0) = 1 — a cos 2 
so that the transformed anisotropy would read 



(2.20) 



t($) = 1 



ma cos 



i - t/2) 



1 + (m — 1) cos 2 



(2.21) 
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III. NUMERICAL INTEGRATION 



We now integrate the rescaled problem, namely Laplace equation for the stream function 
with the boundary conditions Eqs. ( [2.10|) and (|2.16|) . In principle, given an initial condition, 
we should rescale it, evolve it using the rescaled dynamics, and translate the resulting 
interface back to the original variables, but we will not perform any rescaling, since the 
initial condition is free, and the tip-splitting or side-branching character of the result is 
unaffected by the final translation into the original variables. Instead, we will consider the 
rescaled problem in its own, and simulate it by means of the following phase-field model, 



an 

e 2U — = f(9) + e 2 V 2 9 + e 2 K(9)\V9\ + e 2 z ■ (V^ x W), (3.2) 

where 9 is the phase field, e, e are model parameters which must be small to recover the 
sharp-interface equations of the rescaled problem, and we have dropped the tildes of the 
rescaled variables. We have defined f(9) = 9(1 - 9 2 ), and 2^1 = s(9) ■ [V B(6)k{6) + y~ 

k{9) = —V • r(9), with f{9) = and s(9) = f{9) x z. This model was introduced for 



isotropic viscous fingering in Ref. |14| and extensively tested in Ref. ||15|| . From this work 
we know that it will yield converged results for the steady fingers (and, in particular, for 
their widths) if both e < 0.2^/Bo and e < 0.2(1 — c). The only change to be made for the 
anisotropic case is to set B(9) not merely equal to a constant, but to that given by Eq. (|2.17|) 
taking <p = (f)(9) = arccosx • r(9). This gives B(9) = B(<p) + 0(e 3 ), which not only satisfies 
the desired sharp-interface limit, but also ensures that the introduction of anisotropy will 
not result in any extra first-order correction to the free boundary problem, so that the above 
conditions on e, e still hold. 

The same phase-field equations could be used for the circular geometry reinterpreting 
the parameters, since an analogue rescaling yields formally the same result. However, the 
boundary conditions would change. For instance, injection at the center of the cell should 
also be considered. Anyhow, we choose to simulate the well controlled situation in which 
an (unstable) steady finger propagates in a channel of width W. This representation is of 
course exact for single fingers in experiments carried out in the linear Hele-Shaw cell ||, 
but only a (good) approximation for the vicinity of a finger in a multifinger configuration 
(with many fingers) in the circular geometry |7| . 

We investigate the transition between the tip-splitting and side-branching behaviors as 
both the dimensionless surface tension Bo and its anisotropy m — 1 are varied. We are 
interested mainly in the effect of the anisotropy coming from the viscosity (m — 1), so we 
drop that in the original surface tension (a(<p) = a ). The runs use equal viscosities in both 
fluids, c = 0, for reasons of numerical efficiency (e = 0.2), but we do not expect the viscosity 
contrast to affect the stability of the tip for similar reasons for which it does not play any 
role in the (linear) stability of a flat interface. To check this conjecture we ran simulations 
with Bq = 10 -3 both for c = and c = 0.8, two values of the viscosity contrast for which 
a dramatic change in the competition dynamics was seen using the same phase-field model 
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15], and we found that the transition lies in both cases between m = 2 and m = 2.25. These 
c = 0.8 runs are indeed the ones shown in Figs. 1 and 2. 

We use e = 4 x 10 -3 , so that we can simulate accurately with values of the dimensionless 
surface tension down to Bq = 4 x 10 -4 . We first run a steady finger with a large, isotropic di- 
mensionless surface tension B(<j)) = Bq = 10~ 2 , large enough for the finger not to destabilize 
for the amount of (numerical) noise we have and thus reach its steady width and velocity. 
Once this is achieved — see inner interface in Figs. 1(a) and 1(b) — , we perform a "quench" 
in surface tension, i.e., we instantly reduce it to some lower value. Simultaneously, we also 
introduce some amount of anisotropy m — 1. 

The subsequent interface evolution for B = 10~ 3 (and c = 0.8, e = 0.08) is also shown 
within the reference frame moving with the mean interface in Figs. 1(a) (m = 2) and 1(b) 
(m = 2.25) in the form of snapshots at time intervals 0.11. (Simulations used only half of 
the channel and reflecting boundary conditions at its center x = 0). The corresponding y 
position of the interface at the center of the channel (also in the frame of the mean interface) 
is plotted against time in Fig. 2. 

For this value of Bq the finger clearly destabilizes: First its tip widens and flattens (see 
Fig. 1) and therefore slows down (see Fig. 2) for any value of the anisotropy. (Note that 
for t < the tip position would be a straight line in time, since the finger was steady, and, 
in particular, its velocity). Then, for m = 2 the tip continues to flatten and slow down until 
its curvature [Fig. 1(a)] and eventually its velocity in the frame of the mean interface (lower 
curve in Fig. 2) reverse their signs. Finally, the velocity of the interface at the center of 
the channel seems to reach some negative constant value (again, in the frame moving with 
the mean interface) corresponding to the growth of two parallel fingers at each side. We 
identify this reversal of the curvature sign at the center of the finger and this always convex 
tip position vs. time plot with the tip-splitting morphology. 

In contrast, for m = 2.25 the reversal of the curvature sign takes place at some distance 
of the center of the channel, while at the center the curvature increases again [Fig. 1(b)] and 
makes it possible for the tip to speed up again as well, giving rise to a change of concavity 
in the tip position vs. time plot (upper curve in Fig. 2). We identify this reversal of the 
curvature sign at a distance from the center of the channel and this change of concavity in 
the tip position vs. time plot with the side-branching morphology. 

In this way we systematically explore values of the dimensionless surface tension Bq 
ranging from B = 10~ 2 down to B = 4 x 10 -4 . For each value of B we simulate with 
several values of the anisotropy m — 1, and we find that there is a relatively sharp transition 
between the tip-splitting and side-branching morphologies. In Fig. 3 we show for each value 
of Bq (x axis) the two closest values of m— 1 (y axis) for which the two different morphologies 
are observed, namely tip-splitting (circles) and side-branching (triangles). Thus we know 
that the transition line must lie somewhere between the circles and the triangles, and that 
above (larger values of m — 1) and left (lower values of Bq) of that transition line the 
morphology is side-branching, and below and right of it, tip-splitting. This means that the 
critical anisotropy m — 1 above which side-branching replaces tip-splitting decreases with 
decreasing dimensionless surface tension Bq. 

In fact, this critical anisotropy vanishes at Bq ~ 5 x 10~ 4 , and below this value only 
side-branching is observed, even if one uses negative anisotropies down to m — 1 = —0.9, 
which correspond to a viscosity larger in the direction of growth of the finger than in the 
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perpendicular one, and which is not the case of the liquid crystal experiments that motivated 
this study, (m — 1 > — 1 to keep the two viscosities and therefore B((p = 0) finite and 
positive). Of course, the specific value of B for which the critical anisotropy vanishes could 
be affected by the fact that a residual (four-fold) grid anisotropy remains, but it seems 
unavoidable that there is such a (finite) value of Bq, since the transition line curves down as 
Bq is decreased, and for large enough values of B or the anisotropy m—1 the grid spacing 
Ax = e = 4 x is far too fine to affect the effective anisotropy. 

On the other hand, for B > 1.4 x 10~ 3 and for the time elapsed in our runs no clear 
side-branching is actually observed above the transition line extrapolated from lower values 
of i?o, whereas tip-splitting still occurs below that line. For even larger values of B , namely 
Bq ~ 2 x 10~ 3 , not even tip-splitting is observed again within the time elapsed, although 
the steady finger still destabilizes through the widening and flattening of its tip. Finally, for 
B > 10~ 2 the finger is completely stable for the amount of noise we have, as was pointed 
out before. 

IV. DISCUSSION AND CONCLUSIONS 

We have shown that viscous fingering with two different viscosities in two perpendicular 
directions maps to the standard viscous fingering equations (i.e., for isotropic viscosity) with 
an extra two- fold anisotropy in the surface tension, which is such that, together with the 
hypothesis of flow alignment of the director, leads to growth in the direction of maximal 
surface tension. We have simulated the resulting problem using a previously developed 
phase-field model p4]JI3| , and we have found that there is a transition from the tip-splitting 
to the side-branching morphology as either the anisotropy in the surface tension is increased 
or the dimensionless surface tension is decreased. We now draw the connection with the 
liquid crystal experiments of Ref. 0]. 

The observed anisotropy dependence is consistent with the experimental finding that 
there is a transition from tip-splitting to side-branching and back to tip-splitting with tem- 
perature in the nematic phase, since close to the other phases the director alignment, and 
consequently the anisotropy, weaken J7|. The transition is found to be also reentrant with 
injection pressure |7|, which is explained there with the hypothesis that too low pressures 
do not achieve flow alignment, whereas too large ones break down the Hele-Shaw approx- 
imation because of the importance of inertial terms in the hydrodynamic equations, which 
then destroy the flow alignment again. This anisotropy dependence is also consistent with 
simulations of the boundary layer model and the full solidification problem |Hj] , as well 
as with analytical approaches to solidification |§. 

As for the dimensionless surface tension B dependence, one first needs to relate the 
values of B used in the channel simulations to the experimental parameters in the circular 
geometry. To do this we consider a virtual channel whose walls are placed at half the distance 
between a finger and its nearest neighbors. The channel width W is given by this distance 
between adjacent finger tips, whereas the effective injection velocity Voo turns out to be 
the ratio between the injection pressure and R, the mean distance between a tip and the 
injection point. Then, the following dynamic picture of a typical experiment in the circular 
cell emerges: Initially some fingers develop. If the anisotropy, m — 1, is strong enough, their 
tips are stable (which corresponds to a point above the transition line in Fig. 3, where the 
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observed morphology is side-branching). As these fingers grow radially, W increases as R, 
whereas the effective driving force, (771 — 772)^, decreases as 1/R, so that the dimensionless 
surface tension B they experience is found to decrease as 1/R. Thus, the corresponding 
point in Fig. 3 moves to the left (lower values of Bq), and the side-branching behavior is 
preserved. 

In contrast, if the anisotropy m — 1 is not strong enough, the tips split (the corresponding 
point is below the transition line in Fig. 3, where the observed morphology is tip-splitting). 
As a result the number of fingers increases, which then compensates for the growth of the 
distance between finger tips as R in such a way that the effective dimensionless surface 
tension B keeps roughly steady during the pattern development, so that the corresponding 
point in Fig. 3 basically does not move. Thus the transition line is not crossed and the 
tip-splitting behavior is also preserved. In this case, we can estimate B in the experiments 
to be of order 10~ 3 , whereas m is said to be around 2 at the transition 0. This is indeed 
very close to our transition point B = 10~ 3 ,2.05 < m < 2.1 in Fig. 3, but it should be 
taken into account that the value of Bo below which fingers destabilize is known to depend 
on the amount of noise present |]T2[. Accordingly, one expects the whole transition curve to 
be shifted to different values of B Q for a different amount of noise. The latter, however, is 
uncontrolled both in the experiments and the simulations. 

On the other hand, the dependence of the value of the anisotropy at the transition on 
the driving force has been seen in viscous fingering with an etched cell @]. However, we find 
the critical amount of anisotropy for side-branching to replace tip-splitting to vanish at a 
finite value of Bq. Below this value we only observe side-branching. 

A more realistic model of viscous fingering in a liquid crystal should include a velocity 
dependent viscosity. In general the resulting nonlaplacian character of the problem could 
not be avoided, but in principle it would still be possible to simulate the dynamics by means 
of a phase-field model. 
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FIGURE CAPTIONS 



Fig. 1 Destabilization of the tip of a (stationary) finger after instantly decreasing Bq 
to Bq = 1CT 3 at the time of the first interface shown. Successive interfaces are shown in the 
reference frame moving with the mean interface at time intervals 0.11, for c = 0.8, e = 0.08. 
The latest interface is represented in bold, (a) Tip-splitting for m = 2. (b) Side-branching 
for m = 2.25. 

Fig. 2 y coordinate of the interface at the center of the channel (x=0) in the reference 
frame of the mean interface as a function of time corresponding to Figs. 1 (a) (lower curve) 
and 1(b) (upper curve), t = 0(0.88) corresponds to the first (last) interface shown there. 

Fig. 3 Transition between tip-splitting (circles) and side-branching (triangles) as a 
function of the surface tension anisotropy m — 1 and the dimensionless surface tension B . 
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